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ABSTRACT 
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C^I ' Evolution of power spectrum is studied for non-Gaussian models of structure 

• formation. We generalize the dark-matter-approach to these models and find 

pL| ' that the evolved spectrum at weakly nonlinear regime is mainly determined by 

CSl ' a simple integral of the dark-matter-halo mass function in this approach. We 

also study the change of the nonlinear spectrum within the current observational 
^ constraint of the primordial non-Gaussianity. 

^ ■ Subject headings: cosmology: theory — large-scale structure of universe 
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1. INTRODUCTION 

The power spectrum (or equivalently, the two-point correlation function as its Fourier 
transform) is the most fundamental measure to quantify matter clustering. Observational 
■ determination of it will bring us basic information of our universe. For example we can 
make constraints on cosmological parameters or on statistical aspects of the initial matter 



o 



fluctuations that are the seeds of structure formation. Evolution of power spectrum by 
gravitational instability is also an interesting problem of nonlinear physics. 

It is often assumed that the initial density fluctuations obey random Gaussian distri- 
butions. Though this is the simplest assumption from statistical point of views, its origin is 
explained by the simplest inflation scenario. Many analyses of the large-scale structure (LSS) 
in the universe are performed under this assumption. But other models (including some infla- 
tionary models and defect models) predict non-Gaussian initial fluctuations (Vilenkin 1985, 
Vachaspati 1986, Srednicki 1993, Peebles 1997, Linde & Mukahanov 1997) and possibilities 
of these models should be observationally investigated and have recently called much atten- 
tions {e.g. Peebles 1999a, 1999b, Ferreira, Gorski & Magueijo 1998, Robinson, Gawiser & 
Silk 1998, 2000, Koyama, Soda & Taruya 1999, Willick 2000). These recent analyses involve 
the large-scale structure, properties {e.g. time evolution, spatial clustering) of clusters and 
CMB anisotropics. It now becomes clear that defect models cannot generate the observed 
CMB data by their own account but inflation models can explain them well {e.g. Albrecht 



- 2 - 



2000). Recent CMB data by BOOMERanG (Lange et al. 2000) and MAXIMA (Balbi et 

al. 2000) seem to suggest a weaker secondary peak than expected in a simple inflation sce- 
nario. Bouchet et al. (2000) discussed that this might be resolved by a hybrid scenario 
(inflation+cosmic string) in which some component of perturbations come from topological 
defects. 

As for the quantitative measurement of the primordial non-Gaussianity, Contaldi et al. 
(2000) studied the third-order moment of the temperature anisotropics using 4yr COBE 
Differential Microwave Radiometer (DMR) data and found their most conservative estimate 
as {{AT)^) = —6.5 ± 8.7[10^/iK^] (95% CL). Even though the magnitude of the error bar 
is large due to the limited sensitivity and angular resolution, traditional Gaussian model 
(((AT)^) = 0) is consistent with their result. Feldman et al. (2000) analyzed the bispectrum 
of the IRAS PSCz catalog (Saunders et al. 2000) and obtained a constraint for the primordial 
dimensionless skewness as 5" = (5^) / (5^)^^^ < 0.52 (95% GL) for x^-models. Physics of 
GMB are much simpler than that of the galaxy distribution. The signal of the primordial 
non-Gaussianity could be masked by effects of biasing or nonlinear gravitational evolution 
(Verde et al. 2000). Observed data of both CMB and LSS will be significantly improved 
from MAP and Plank satellites, SDSS and 2dF surveys. Verde et al. (2000) insisted that 
GMB would be a better probe of the primordial non-Gaussianity for several non-Gaussian 
models and LSS data would be useful to constrain the biasing if primordial Gaussianity is 
suggested by GMB map. 

So far nonlinear effects of matter clustering from non-Gaussian initial conditions are 
studies by numerical simulations (White 1999, Robinson & Baker 2000 and references therein) 
or by perturbative methods {e.g. Fry & Scherrer 1994, Ghodorowski & Bouchet 1996, Scoc- 
cimarro 2000, Durrer et al. 2000). In this article we study evolution of the matter power 
spectrum for non-Gaussian models using the dark-matter-halo approach that can be traced 
back to Peebles (1974) and McGlelland and Silk (1977) (see also Scherrer & Bertschinger 
1991). For Gaussian models this approach is confirmed to reproduce numerical results very 
well from linear to nonlinear scales and regarded as an excellent tool to explore the matter 
clustering. It seems interesting to see how the nonlinear spectrum can be different within 
the currently allowed region of the primordial non-Gaussianity. 

Many works have been recently performed with this approach for Gaussian models. 
For example nonlinear evolution of the power spectrum or the bispectrum of density field 
is studies by Ma & Fry (2000a, 2000b, 2000c), Seljak (2000), White (2000), Peacock & 
Smith (2000), Scoccimarro et al. (2000) and weak lensing effects by Cooray, Hu & Miralda- 
Escude (2000) and Gooray Sz Hu (2000). Komatsu & Kitayama (1999) investigated the 
Sunyaev & Zeldovich effect using a similar method. Here we simply extend this approach 
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to non-Gaussian models. This approach is made from three basic ingredients (i) the mass 
function of dark-matter-halos, (ii) the bias parameter of halos relative to the large-scale 
density fluctuations and (iii) the density profile of halos. Numerical analysis (Robinson 
& Baker 2000) supports that a straightforward extension of the Press & Schechter (1974) 
mass function is effective for non-Gaussian models. But validity of simple extensions of two 
other elements have not been checked numerically. Therefore we make a careful analysis and 
limit our investigation in the range where results are expected not to depend largely on the 
details of these two elements. Our goals in this article are set to the following two points. 
Firstly, we intend to roughly understand the effects of primordial non-Gaussianity on the 
weakly nonlinear evolution of the power spectrum. Secondary, we estimate the difference 
of the evolved spectra within the current observational constraint of the primordial non- 
Gaussianity. 

This article is organized as follows. In §2 we describe basic properties of the dark- 
matter-halo approach with its basic three ingredients, namely, the mass function of halos 
(§2.1), the bias parameter (§2.2) and the density profiles of halos (§2.3). In §3 we analyze toy 
models to clarify the general effects of non-Gaussianity on the weakly nonlinear evolution of 
power spectrum. In §4 we investigate realistic models. We study a flat A-dominated model 
{flo = 0.3, Ao = 0.7 and h = 0.7) with a fixed linear CDM spectrum and evaluate evolution of 
power spectra for log-normal probability distribution functions that are motivated by work 
of Robinson & Baker (2000). §5 is devoted to a brief summary. In appendix A we discuss 
the scale dependence of the skewness parameter for x^-models. 



2. DARK MATTER HALO APPROACH 

In the dark-matter-halo approach the nonlinear density field is given by superposition 
of halos with various masses. The density profile of each halo is assumed to be determined 
by its mass. The nonlinear power spectrum PNiik) is constituted by two terms as 

PNL{k)^P2h{k)+Pih{k). (1) 

The two-halo term P2h{k) counts contributions of two points coming from two different halos 
and is given by 

2 

PL{kl (2) 

where Piik) is the linear power spectrum. The one-halo term Pih{k) counts two particles 
within the same halo 

Pih{k) = p^^ / dmN{m)m^u{k,mf. (3) 



P2h{k) = 



dmN{m)mb{m)u{k, m) 
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In the above equations pb is the background density of the universe and N{m)dm is the 
number density (mass function) of halos with mass from m to m + dm. The factor h{m) 
is the bias parameter and represents distribution of halos relative to the large-scale density 
fluctuations. The function u{k,m) is the Fourier transform of the density profile of a dark- 
matter-halo with mass m. The original form of the dark-matter halo approach was proposed 
by Peebles (1974). This approach has recently called much attention as our understandings 
of its ingredients have largely developed. In the following three subsections we study these 
basic ingredients. 



2.1. Mass Function 

Wc use the Press & Schcchtcr (1974) formalism for the mass function N{m) of collapsed 
objects. This formalism was first given for Gaussian random fields but can be straight- 
forwardly extended to general models {e.g. Lucchin & Matarrese 1989, Chiu, Ostriker & 
Strauss 1998, Koyama, Soda & Taruya 1999, Robinson, Gawsier & Silk 2000, Willick 2000, 
Matarrese, Verde & Jimenez 2000). For simphcities we assume that the hnear one-point 
probability distribution function (hereafter PDF) p(5, a) for the smoothed density field is 
written as 

p{5,a)d5^a-^p{S/a)d5, (4) 

which means that the shape of the PDF is scale invariant. This ansatz is correct for the 
random Gaussian fluctuations, but not guaranteed for general non-Gaussian models. We 
comment on this scale dependence in §4 and appendix A. 

The volume fraction of points that belong to collapsed objects with mass larger than m 
is given as follows 

poo 

/ p{u)du, (5) 

•J 5c/cT{m) 

where 6c — 3(127r)^/^/20 = 1.69 is the critical linear density contrast for the spherical top-hat 
collapse and a{m) is the root-man-square mass fluctuation in a sphere of mass scale m. F 
is the normalization factor to account for all masses in the universe 

POO 

F= / p{iy)diy. (6) 
Jo 

For a symmetric profile p{-~i') = p{i^) we have F = 1/2 (see also Bond et al. 1991, Lacey & 
Cole 1993). From equation (5) the number density of halos with mass from m to m -|- dm is 
given as 

1 -1 *9 1 _i d\na(m) 
N[m)dm = —F p^m — — / piujdu = —F pbpii'jm v ; , (7) 

J5,/a{m) dm 
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where we have denoted the normahzed density contrast v = 5c/ a{m) . We define the nonhnear 
mass m^: by equation (T(m*) = 5c- For power- law models with Piik) oc A;" the linear mass 
fluctuation a{'m) is given as 

a{m)^5cl — ] . (8) 

In this subsection we have followed the basic analysis of the Press & Schechter method. 
For random Gaussian fluctuations some refinements have been proposed to reproduce N- 
body data better {e.g. Sheth & Tormen 1999, Jenkins et al. 2000). If we write down the 
mass function in the following form 

iV(m) = -A\[^^^ue-'^-''/\l + {au'm- (9) 
V TT ammm^ 

The original Press & Schechter function corresponds to a) = (0.5,0,1). The mass 

function given by Sheth & Tormen is given by parameters {A,p, a) — (0.322, 0.3, 0.707). 

Robinson & Baker (2000) have recently confirmed that the simple extension of the Press 
& Schecheter formula (7) can accurately fit the mass function of clusters for several non- 
Gaussian models with typical error ~ 25 percent (see also Avehno & Viana 2000). This 
result encourages our analytical study. 



2.2. Bias Pctrameter 

Next we study the bias parameter b{m) using (i) the peak-background splitting method 
(Kaiser 1984) and (ii) a simple extension of the Press & Schechter method (Mo & White 
1997, Mo, Jing & White 1997, Sheth & Lemson 1999, Sheth & Tormen 1999). When a large- 
scale perturbation 5ls is added to the local (small scale) density field (^55 (a?), the density 
contrast for 5ss{x) that is critical to the spherical top-hat collapse becomes 6c — Sls- Then 
the Eulerian bias parameter for halos with mass m becomes (Mo & White 1997, Robinson, 
Gawiser & Silk 1998, Koyama, Soda & Taruya 1999) 

, d\nN(m) 1 fp'M \ , , 

From equations (7) and (10) we can easily confirm the following equation ^ 

t-oo 

p^^ I b{m)N(m)mdm — 1. (11) 







^This equation holds for general non-scaling PDF in the form p{6, a). 



-6- 



As we see in §4 this equation is useful to know the behavior of the nonhnear power spectrum 
-FWl(^) at large scale (small k). There are several analytical works about the bias parameter 
for non-Gaussian models (Robinson, Gawiser & Silk 1998, Koyama, Soda & Taruya 1999). 
But validity of the formulas (10) for a non-Gaussian PDF have not been checked with 
numerical simulations. But we expect that it works well, as this formula is simply derived 
from the Press & Schechter method that has been checked numerically for some non-Gaussian 
models as commented in the previous subsection. 

For the mass function (9) given for random Gaussian fluctuations we have 

, / \ — 1 2p ^, 



2.3. Density Profile of Dark Matter Halo 

Recent numerical simulations suggest that the average density profile of cold-dark-mater 
(CDM) halos has nearly universal shape pm{i^) characterized by its mass m (Navarro, Frenk 
& White 1996, 1997, Moore et al. 1998, 1999, Kravsov et al. 1998, Fukushige & Makino 
1998, Jing & Suto 2000). These numerical simulations are performed with Gaussian initial 
conditions, but here we assume that the density profile depends only weakly on the linear 
PDF as the halo profile would strongly reflect nonlinear gravitational dynamics rather than 
the initial conditions. This assumption may be a potential weakness of our dark-matter- 
halo approach especially at strongly nonlinear regimes. However as we see in §4, the power 
spectrum PNhik) at weakly nonhnear regimes does not depend details of the halo profile. 
In this article we adopt the universal profile of dark matter halo given in Navarro, Frenk & 
White (1997) as follows 

/ \ Pmy) /ION 

Um[r) = (13) 

m 

^ (14) 



vir I 



where R^r is the virial radius of the halo. The halo mass m is written by Rmr as m = 
47ri?^j^PbA/3 with the average density contrast A = 200 (Einstein de- Sitter model) and 
A = 340 (fiat model with VLq = 0.3 and Aq = 0.7). The concentration parameter c determines 
the transition of inner and outer regions of a halo and depends on its mass m, roughly 
reflecting the formation epoch of the halo (Ma & Fry 2000b). We have c = (3{n)a{m) with a 
coefficient (5{n) that depends on the linear power spectrum. The normalization factor / in 
equation (14) is given by 

/=(ln(l + c)-c/(l + c))-^ (15) 
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The density profile (14) in central region becomes p{r) oc r°'. The original result in Navarro, 
Prenk & White (1996) was a = —1. But some other simulations show different results. For 
example, Moore et al. (1998, 1999) found a = —1.5 and Jing & Suto (2000) argued that 
the index a is not a universal value but depends on the mass scale of halos. As we see later, 
the nonlinear evolution of power spectrum in weakly nonlinear regime is not sensitive to the 
details of the halo profile (see also Seljak 2000). Thus we use a = —1 in this article. 

Finally the Fourier transformed function u{k, m) is defined in terms of the real space 
density profile Umir) as 

u{k,m) = j dp X GXY>{—ix ■ k)ujn{r) (16) 

/Sill /u?^ 
—g—r'^Um{r)dr. (17) 

For the profile (14) an explicit analytic formula of u{k,m) was given by Scoccimarro et al. 
(2000). At small wave number k the above Fourier transformation becomes a simple volume 
integral and we have u{0,m) = 1 due to the definition of Um{r) given in equation (13). 



3. TOY MODELS 

In this section we discuss effects of initial non-Gaussianity on the evolution of the power 
spectrum in Einstein de-Sitter background. We intend to extract out their general relations. 
From this standpoint we do not try to mimic realistic PDFs so seriously but study toy models 
that arc different from Gaussian distribution. We also limit our analysis to scale-free initial 
spectra. With results obtained in this section we study models relevant for observational 
cosmology in §4. 



3.1. Density Distribution Functions 

We use the following four PDFs for the linear smoothed density fields. As commented 
in §2.1 we assume that the PDFs are scale invariant. The first example is the standard 
Gaussian model (G-model), 

Pg{'^) = A— ^-^p(~'Tr] — oo < < oo. (18) 



V2^ "V 2 

This PDF is symmetric around the mean density u — and we have F — 0.5. Here the 
factor F is defined in equation (6). 
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The second example is the Chi-square (x^)-model (C-model) with two degrees of free- 
dom. A x^-model is often adopted as a typical non-Gaussian model and its origin in cos- 
mological context was recently proposed by Peebles (1999a) using an inflation model. The 
PDF for the unsmoothed field is 



exp(— i/ — 1) (—1 < u < oo), 
{u < -1), 



(19) 



with F — 1/e. We should notice that the PDF for the smoothed density field is same as the 
above unsmoothed one. This point is discussed in the next section and appendix A. 



Next example is a log-normal distribution. This distribution is characterized by a pa- 
rameter a{> 0) as follows 

<|ln(Vl + Q;^(l + i^tt))|^ 

: CXp 



a 



PLiv(i^, a) = < (l + z/a)^27rln(l + a2) 



21n(l + a2) 







—a ^ < u < oo, 
(20) 

At a = wc recover the Gaussian distribution pciy) = Pln{j^, 0). The normalization factor 
F for this distribution is given in terms of the error function as 



-erfc 
2 



1 Jln(l + Q;2) 



(21) 



where the error function is defined as erfc(a;) = ^/2/7^ exp(—y'^/2)dy. 

Here we add two models generated from equation (20). We define the LNl-model by 



PLNli^) = Pln{^, 1) 



(22) 



This profile has larger positive tail than the G-model. We have F = 0.3386. The final 
example (LN2-model) is obtained by a simple replacement of the previous model as 



PLN2{i^) =PLNli-J^) ^PLNi-T^, 1)- 



(23) 



Note that we have PLN2W) = for i/ > 1. Existence of halos with mass larger than m* is 
prohibited in this model. As we see in the following section, the shape of the positive tail 
is very important to study nonlinear evolution of the power spectrum. The probability for 
z/ > 3 is 1.34 X 10-3 (G-model), 1.83 x lO'^ (C-model), 1.87 x lO'^ (LNl-modcl) and 
(LN2-model). For z/ > 5 we have 2.87 x 10"^ (G-modcl), 2.48 x lO"-'^ (C-model), 5.11 x lO^^ 
(LNl-model) and (LN2-model). We can easily confirm that all of these models have 
vanishing means 

(i/) = / up{u)dv = 0, (24) 
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and are normalized properly 

/oo /»oo 
p{u)du = 1, = / u^p{u)du = 1. (25) 

oo J —oo 

The skwness parameter 5" = (i/^) / (i/^)^^^ is a fundamental measure to characterize 
PDFs. 2 We have S^O (G-model), 2 (C-model), 4 (LNl-model) and -4 (LN2-model). For 
the general log-normal distribution Pln{j^, we have 

SLN = a{3 + a^). (26) 

We use this relation in section 4. 



3.2. Results 

Using ingredients described in §2 and various linear PDFs given in §3 we calculate the 
nonlinear power spectrum PNhik) with equations (1), (2) and (3). We use power- law models 
for initial spectra 

PL{k)^Ak^, (27) 

where A is the normalization factor but irrelevant for present analysis. Time evolution of 
power spectrum can be reduced to rescaling of spatial length due to our scale-free initial 
conditions. We introduce the nonlinear wave-number kj^L that is defined by 



1 rkNL 

J PL{k)k^dk - 1. (28) 



This wave-number is the most fundamental scale for the present analysis. In this article we 
call the weakly nonlinear scale for the wave number k 

0<k< kNL- (29) 

Wc first calculate the nonlinear power spectra PNiik) for random Gaussian fields and 
compare our analytical predictions with fitting formula that is obtained from N-body simu- 
lations with the same Gaussian initial conditions (Peacock & Dodds 1996, see also Hamilton 
et al. 1991, Jain, Mo & White 1995). We fix the coefficient (3{n) ~ 3 for the concentra- 
tion parameter c(m), but our results are insensitive to this parameter. It is found that our 



^Definition of skewness parameter is different from that usually used for Gaussian initial condition 
(Bouchet et al. 1992). 
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model cannot reproduce the numerical results for spectral index with n = in contrast to 
the case with n < — 1. We can see this tendency from figures l.a to I.e. It is not clear 
why this dependence appears. Note that the evolved spectrum (thick solid line) for n = — 1 
obtained from the fitting formula is smaller than the linear spectrum (dash-dotted line) as 
predicted by the second-order perturbation theory {e.g. Makino, Sasaki & Suto 1992). In 
the framework of the dark-matter-halo approach this fact might be related to the possibility 
that the total halo mass fraction does not converge to unity (Jenkins et al. 2000) and the 
constraint (11) does not hold. As we see below, without the constraint (11) the two-halo 
term P2h{k) can be smaller than the linear spectrum Piik) with nonnegligible amount at 
weakly nonlinear scale. 

Here we limit our analysis for spectral indexes with n < —1 where our approach is 
expected to work well. These are close to the effective slope of a typical CDM power spectrum 
at weakly nonlinear regime c < 1. In figures l.a to l.c we plot the nonlinear power 
spectrum PNiik) in the form A(/c) defined by 

A(*) . (30) 

for spectral indexes n = —2, —1.5 and —1. In these figures contributions from the one-halo 
term Pih{k) and the two-halo term P2h{k) are presented separately. Note that our analytical 
prediction for the Gaussian models reproduce the results from N-body simulations well in the 
range k < k^L for n < —1.5 This fact indicates validity of the dark-matter-halo approach. 
For n = —2 model we also evaluate the nonlinear spectrum with using generalized formulas 
(9) (12) as in Seljak (2000). Our results shows good agreement with this results, in figures 
l.a to l.c w have used the fitting formula of Peacock & Dodds (1996). We find that the 
fitting formula of Jain et al. (1995) agrees with our analytical result better for n = —2 model 
but worse for n = —1 model. 

Next wc discuss behavior of the nonlinear power spectrum PNiik) for various PDFs 
given in §3.1. We present numerical results in figures l.a to l.c. These are obtained by 
fully using the dark-matter-halo approach. In this calculation we evaluate all the three 
elements, the mass function N{M), the bias parameter b{M) and the halo profile u{k,M). 
But information of the mass function is most important at weakly nonlinear regime, as we 
see below. 

For all models the two-halo term in the range k < /catl is almost same as the linear 
power spectrum Piik), as shown on figures for the two-halo term. In this regime we have 
u{k, m) ~ 1 in the dominant contribution of the integral in equation (2) (see also figure 2). 
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Then we obtain the following relation from equation (11) 

/*oo /»oo 

/ dmN[m)mh{rn)u{k,rn) c::^ / dmN{m)mb{m) — pb, (31) 
Jo Jo 

and 

P2h{k)^PL{k), (32) 
(see also Sheth et al. (2000) for a same kind of analysis for the two-point correlation function). 

Therefore we can roughly understand the weakly nonlinear evolution of power spectrum 
only by analyzing the one-halo term 

Pih{k) = Pb^ J m^N{m)u{k, mfdm. (33) 

As shown in figures l.a to l.c, the contribution of A(A;) from one-halo term Pih{k)k^ / (27r^) 
is nearly on a straight line proportional to k^. This means Pih{k) — const and indicates that 
/c-dependence caused by the halo profile u{k, m) is weak. Thus we approximate the one-halo 
term at weakly nonlinear regimes k < k^i by setting u{k, m) — 1 

Pih{k) ^ Pih{Q) = pf J dmN{m)m\ (34) 

We have a constraint for the function mN{m) due to the conservation of mass as 

dmN{m)m — pf,, (35) 



/ 

Jo 



and Pih{0) in equation(34) becomes larger for a mass function N{m) with more abundant 
massive objects. Thus we can roughly understand the weakly nonlinear correction of power 
spectrum using information of mass function in the form given in Pih{0). 

Note that the present analysis at weakly nonlinear regime k < k^L based on equations 
(32) and (33) does not depend on details of the halo density profile. Therefore in this regime 
we should not be troubled seriously with potential problems about the reliability of the uni- 
versal profile for non- Gaussian models. This also means that our analysis is not sensitive 
to the bias parameter b{m) (cq.[10]). In the region where u{k,m) = 1 is a good approxi- 
mation, we need the bias parameter only in the form of equation (11). But this equation 
represent a simple requirement that the large-scale fluctuation is analyzed perturbatively 
PNL{k) ~ PL{k). Here we should notice that we have not taken into account the nonlinear 
biasing effects or the exclusion effects of halo-halo clustering. These effects might be im- 
portant in weakly nonlinear regime (Scoccimarro et al. 2000), but we neglect them in this 
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article. At present there is no simple prescription to deal with these two opposite effects. 
This is a generic problem of the dark-matter-halo approach. 

Using results for n = —2 model we discuss effects of the linear PDF on the weakly 
nonlinear evolution of the power spectrum. As shown in figure l.a the nonlinear correction 
for the LN2-model is very weak. In this model there are no halos with m > m* (see figure 
2) and the integral / N{m)m'^dm is much smaller than other models. The weakly nonlinear 
effect is strongest for the LNl-modcl as expected from its mass function in figure 2. Even 
at k k^L magnitude of the nonlinear power spectra can differ by a factor of 10 depending 
on the initial PDFs. 

Next we compare the behavior of nonlinear power spectrum PNiik) for different spectral 
indexes n. As n becomes larger, the nonlinear correction becomes smaller (figures la-lc). 
All nonlinear power spectra P]\fL{k) for n = — 1 model are very close to the linear spectrum 
Piik), but these for n = —2 model depend largely on the adopted PDF. This n-dependence 
can be also understood by using the integral Pi/i(0) given in equation (34). We have the 
following simple relation for the root-mean-square mass fluctuations a{m) 



As the index n increases, abundance of large mass objects becomes smaller. Then the 
magnitude of the one-halo term ~ Pi/i(0) decreases and we have smaller nonlinear corrections. 



In this section we study nonlinear evolution of the power spectrum for realistic models. 
We fix the cosmological parameters at h = 0.7, Qq = 0.3 and Aq = 0.7 that are compatible 
with CMB {e.g. Lange et al. 2000) and high-redshift SNe data (Ricss et al. 1998, Perlmutter 
et al. 1999). We use the primordially scale invariant {n = 1) spectrum with CDM transfer 
function from Bardeen et al. (1986). The power spectrum is COBE- normalized and we have 
as — 0.9. For these parameters the nonlinear wavenumber k^L defined in equation (28) 
becomes kNL — 0.19Mpc. 

Robinson & Baker (2000) studied evolution of the cluster abundance in representative 
non-Gaussian models, such as, Peebles isocurvature cold-dark-matter (ICDM) model or cos- 
mic string model. This ICDM model corresponds to a x^-model with one degree of freedom 
(see appendix A). Peebles (1999b) show that the power spectrum for this model with param- 
eters Qo = 0.2, A = 0.8 and = 2.4 (see eq. [A14] for definition of n<^) is in good agreement 
with observed spectrum of CMB and the large-scale structure. 




(36) 



4. 



CDM SPECTRUM 
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Robinson & Baker (2000) also analyzed PDFs of smoothed linear density fields for these 
non-Gaussian models and found that the PDFs are well fitted by the log-normal distributions 
that are nearly scale independent at weakly nonlinear scale. For example, they did not detect 
scale dependence for the ICDM model and the string+HDM model (their table 1). Thus, 
motivated by their work, we investigate the log-normal PDFs given in equation (20) and 
neglect the scale dependence. We characterize the non-Gaussianity of a log-normal PDF by 
the skewness parameter S = (5^) / {6^)'^^'^ (see eq.[26]). For models in Robinson & Baker 
(2000) we can estimate the parameter S from their table 1. We have S ~ 2.2 (somewhat 
smaller than result ~ 2.5 given in Peebles 1999b) for the ICDM model with = —2.4 
and S ~ 0.5 for the string-|-HDM model. The parameter 5" for the ICDM model is larger 
than the observed constraint S < 0.52 (95% CL) given in Feldman et al. (2000). A simple 
improvement is to increase the number of scalar fields. This number corresponds to 
the degree of freedom in x^-model. As discussed in appendix A, the parameter S scales are 
S oc N~^^'^ and the required number of freedom N is N > (2.2/0.52)^ ~ 20 (Feldman et 
al. 2000). 

In figure 3. a we show the nonlinear spectrum at the wave number k = k^L as a function 
of the parameter S. The concentration parameter is fixed at c(m) = 9(m/m*)~°'^^ (Bullock 
et al. 1999). We first calculate the ratio 

PNLjk) ^ P2h{k) + Pihjk) 

Pdk) - PL{k) ' ^ ^ 

in the full dark-matter-halo approach. We also evaluate the simple approximation with 
putting u{k, M) — 1, namely 

PM!^^ P^ (38) 

PL{k) - + PL{k) ^^^^ 

to check the effects of the bias parameter h{M) and the halo inner density profile u{k,M). 
As commented before, validity of these two elements for non-Gaussian models have not been 
clarified numerically so far. The simple approximation (38) shows reasonable agreement with 
the full result (37) and we can expect that our analysis is not seriously affected by uncertainty 
of two elements in the framework of the dark- matter-halo approach. Figure 3. a shows that 
nonlinear spectrum depends strongly on the primordial non-Gaussianity measured by the 
parameter S. For the Gaussian model (-S" = 0) we have Pnl/Pl ~ 1-3 but Pnl/Pl ~ 2 for 
5" ~ 1. For the value S — 0.52 (on 95% CL of Feldman et al. 2000) the correction becomes 
Pnl/Pl < 1-55 and about 20% larger than the Gaussian model. We also calculate the 
ratio Pnl/Pl for the Gaussian model using the fitting formula of Peacock & Dodds (1996). 
The result is ~ 1.2 and close to our result 1.3. 



In figure 3.b we show the nonlinear spectrum at the wavenumber k — 0.66A;jvl where 
the hnear CDM spectrum becomes Pz;,(0.66/cArL) = '^PLik^L)- Apparently the correction 
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becomes smaller than at k — k^L- We can easily understand this fact. At weakly nonlinear 
scale we have 

PNL{k) - PL{k) ~ Pi^O) (39) 

and the nonhnear correction behaves as Pnl/Pl — 1 + -Pi/i(0)/Pl(^) in the dark- matter-halo 
approach. 

5. SUMMARY 

In this article we have studied evolution of the matter power spectrum in non-Gaussian 
models of structure formation. We have used the dark-matter-halo approach that is an 
excellent method for analyzing matter clustering from linear to nonhnear scales (Seljak 2000, 
Ma & Fry 2000b). This approach contains three basic ingredients (i) the mass function N[m) 
of dark-matter-halos, (ii) the bias factor b{m) of halos relative to the large-scale density 
fluctuations and (iii) the density profiles u{k,m) of halos. Recent numerical simulations 
(Robinson & Baker 2000) have suggested that a simple extension of the Press & Schcchter 
formula provides good fits to mass function of clusters also in non-Gaussian models. For the 
bias factor b{m) and the density profiles u{k, m) it has not checked numerically whether we 
can simply extend the formulas of Gaussian models to non-Gaussian models. Therefore we 
take a moderate position to limit our investigation up to weakly nonlinear regimes k < k^L- 
In these regions our analysis supports that the evolved power spectrum does not depend on 
details of the bias factor b{m.) or the density profile u{k, m). 

The nonlinear power spectrum predicted by the dark-matter-halo approach is consti- 
tuted by two terms. The two-halo term counts two points in different halos and the one-halo 
term represents contributions of two points within same halos. We have shown that the 
two-halo term is almost same as the linear power spectrum in weakly nonlinear regimes and 
nonlinear correction mainly comes from the one-halo term. This term is roughly approxi- 
mated by the integral of mass function N{m) as J N{m)m^dm. Using this integral we can 
understand various aspects of the weakly nonlinear effect. For example, nonlinear correction 
is larger than a PDF with more positive tail or smaller spectral index n. 

We also study a COBE-normalizcd CDM spectrum in a flat model {VLq = 0.3, Aq = 0.7 
and h = 0.7) with log-normal PDFs. This shape of PDF is motivated by work of Robinson 
& Baker (2000) who found numerically that PDFs of many representative non-Gaussian 
models are fitted well by log-normal distributions. We show that the nonlinear spectrum 
Pnl a,t k — kNL — 0.2Mpc~^ becomes Pnl/Pl — 1-3 for the Gaussian model but ~ 1.55 
for the log-normal model with S — 0.52 that is on 95% CL of the observed constraint by 
Feldman et al. (2000). These results suggest that the ratios Pl/Pnl for ACDM model at 
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k < kNL — 0.2Mpc ^ could not largely (say factor of 2) vary within currently allowed region 
of the primordial non-Gaussianity. 

The author would like to thank an anonymous referee for valuable comments to improve 
this manuscript. My work is supported by Japanese Grant-in-Aid No. 0001416. 



A. SCALE DEPENDENCE OF THE SKEWNESS IN x^-MODELS 

The scale dependence of a non-Gaussian PDF becomes highly comphcated as it is gen- 
erally very difficult to deal with infinite degrees of freedom of nonlocal quantities. Here we 
study the scale dependence of the PDF using the skewness parameter for Xn model. 

Let us discuss Xn distribution with N degrees of freedom. The unsmoothed density 
field p{x) is written by N independent Gaussian random fields 0i(a;) {i — 1, - ■ ■ ,N) 

N 

p{x)^Y.^,{xf. (Al) 

i=\ 

We assume that these N fields have same statistical character. Moments of the unsmoothed 
field p{x) is given in terms of the variance of the basic Gaussian fields <pi (the bracket 
(•) represents ensemble average). Wc can easily calculate their moments as 

{p{x)) = N{ct>l), {5p{xf)^2N{ct>l)\ 
{5p{xf) = %N{4>1)\ {5p{xY)^{A%N + l2N'){<t>l)\ (A2) 

where we have defined the overdensity 5p{x) hj 5p{x) = p{x) — {p{x)). Furthermore we can 
derive the explicit formula for the PDF of p as 




(1 + z/V2iV-i)^/^-^ 
^^(^)^^ = (2iv-i)(A^-i)/2r(Ar/2) """"P I I 1 + l/ 1 H^' (^3) 

where we have defined the normalized overdensity v hy v = 5p/ {Sp'^)^^'^ and r(-) is the 
Gamma function. As the number N increases, the PDF becomes closer to the Gaussian 
distribution. This seems reasonable considering the central limit theorem. For smaller 
degrees of freedom N, the PDFs show strong non-Gaussianity. 

The skewness parameter 5" is a fundamental quantity to characterize the non-Gaussian 
profile. For x^-model we have the following expression from equation (A2) for moments 
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Note that using the skewness parameter S, we can specify the x^-distribution that is com- 
pletely characterized by the degrees of freedom N. 

The simphcities of relations given so far are mainly due to local nature of the density 
field p{x). Every statistical characters of the field are reduced to simple Gaussian statistics 
of the basic fields (j)i{x). When we discuss smoothing effects on the density field p{x) that 
is an nonlinear combination of Gaussian fields (pi, the situation becomes highly complicated. 
First we define the smoothed density field Pr{x) by 

pji{x) = J dx' p{x')w{x' -X : R), (A5) 

where w{x; R) is a filter function with smoothing radius R. In this appendix we mainly use 
the Gaussian filter 

w[x : R) — — -^^=^= exp 



We expand the fields 4>i{x) in terms of the Fourier modes 4>i{k) 

M^) = j^pj^ J dk(l)i{k) exp(ifc ■ x). (A6) 

As the basic fields 0j obey independent random-Gaussian distributions, their statistical na- 
ture is completely determined by their power spectrum P</,(/c) defined by 

= P4,{k)S,,SD{k + I), (A7) 

where Sij is the Kronecker's delta, is the Dirac's delta function and we have assumed 

that fluctuations are isotropic. Similarly the smoothed density field Pr{x) is written by 4>i{k) 
as 

N 

Pr{x) = ^ j dkdlct>i{k)Ul) exp(i(fc + I) ■ x)W{\k + l\R), (A8) 

where W{kR) is the Fourier transformed filter function of w{x : R). For the Gaussian filter 
we have W{kR) = exp(— /c^i?^/2). With expression (A8) we can write down the second- and 
third-order moments for the smoothed density field in terms of the power spectrum P^{k). 
After some algebra we arrive at 

{5pl) = dkdlP^{k)P^{l)W{\k + l\R)\ (A9) 

{5p\) = J dkdldniP^{k)P^{l)P^{m)W{\k + l\R)W{\l + rn\R)W{\m + kmiO) 

Using these equations we evaluate the skewness parameter Sr for the smoothed density 
field. We can easily confirm that the unsmoothed value S — -^/^/N is recovered by setting 
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R — in these equations. The smoothing effect changes the skewness parameter Sr from 
the unsmoothed value S — ^/WfN. To discuss quantitative effects of smoothing we define a 
parameter defined by 

(Ml) 



Note that the parameter Jr does not depend on the degrees of freedom A^. The effective 
degrees of freedom N' determined form the smoothed skewness parameter Sr becomes 
times original value. 

Next let us evaluate the smoothed skewness Sr numerically and the scaling parameter 
Jr for a sequence of matter fluctuations. We adopt the following from for the power spectrum 
of the Gaussian fields (pi 

P^k) = Ak""^ exp(-A;V2). (A12) 

As the normalization factor A is irrelevant for our linear analysis we simply put A = 1. Due 
to the Gaussian cut-off exp(— A;^r^), small scale power at wave-number k > is strongly 
suppressed and the fields (f)i{x) and p{x) are very smoothed at spatial scale smaller than r. 
This means that the smoothed density field Pr{x) is close to the unsmoothed field p{x) for 
smoothing radius R with R < r. Therefore radius r represents the spatial scale where the 
PDF of the smoothed field pR would become close to the original Xn distribution and the 
scaling factor would be /i? ~ 1. 

Using equation (A9) for variance of the smoothed density field, we obtain the following 
result for the Gaussian filter 



i^Pl) = I ^W^^^exp[-(^l + ^J(.^ + t2) 



sinh(2st). (A13) 



At large smoothing radius R we have {{Sp)"^) oc ^ ^"'^ and the matter power spectrum at 
small wave-number would behave as (Peebles 1999b, White 1999, Scoccimarro 2000) 

Pp{k) oc = (A14) 

where we have defined the spectral index Up of density field by Up = 2n^ -\- 3. 

The expression (AlO) for the third-order moment is a nine-dimcnsional integral. But 
using symmetry with respect to variables fc, I and m, we can simplify the expression. After 
some algebra we arrive at 

jYj^— 9— 3n^ roo roo poo ptt /-tt 

{^Pr} — ;r~6 / '^'^ "^W '^^ / da sin a / db sin b{stu)^~^"''i'Io{tu sin a sin b) 

Jo Jo Jo Jo Jo 

X exp — ( i -h — ) [s^ + t^ + u^) — [St cos a + tu cos a cos b + us cos b ) , (A15) 



— f 1 -|- j (s^ + t'^ + u^) — {st cos a + tu cos a cos b + us cos b) 
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where Io{x) is the 0-th modified Bessel function. Prom equations (A13) and (A15), it is 
apparent that the factor /r is written in the form 

fn = f{R/r). (A16) 

In the large scale limit — > oo, the small scale Gaussian cut-off in the power spectrum 
P<p{k) (eq.[Al2]) would have no effects on the second- and third-order moments as expected 
from equations (A13) and (A15). In this limit we can put r — and our results should 
coincide with results for pure power-law models Pcf,{k) oc A;"-^ (Peebles 1999b, Scoccimarro 
2000). Simple analytic expressions for the matter power spectrum Pp and bispectrum Bp for 
— —2 model are given by Scoccimarro (2000) as follows 

P,(fc) = ^, Bp{k,l,m) = ^-^. (A17) 

With the latter expression we can numerically study the third-order moment {{5pr)^) at 
— > oo. In this case we need only two-dimensional numerical integration that is much 
simpler than the five-dimensional integration (A15). We find the numerical result ^ 

/(oo)=2.6 (n^ = -2). (A18) 

This asymptotic value can be used to check our five- dimensional numerical integration re- 
quired to evaluate the smoothed skewness at an intermediate scale. 

In figure 4 we show scale dependence of the factor f{R/r) for power spectrum (A12) 
with various spectrum indices, = —1, = —1.5, = —2 and = —2.4. These indices 
correspond to Up = 1, Up = 0, Up = —1, and Up = —1.8, respectively. The last choice = 
—2.4 is same as Peebles (1999b). As expected, we have f{R/r) ~ 1 for i? < r reflecting 
the original profllc of the unsmoothed PDF. For a larger smoothing radius R > 20r the 
parameter f{R/r) relaxes to an asymptotic value f{oo) that does not depend on the small 
scale cut-off as explained earlier. The asymptotic value for — —2 model is close to the 
value given in equation (A18) and show the validity of our numerical integration. Peebles 
(1999b) calculated the skewness parameter smoothed by the top-hat filter for a model with 
one degree of freedom = 1 and spectral index at = —2.4, using Monte Carlo integration. 
His result is Sr = 2.46 and corresponds to Jr ~ 1.32. This value is close to our asymptotic 
value /(oo) ~ 1.37 obtained for the Gaussian filter. 

As shown in figure 4, the factor f{R/r) largely depends on the spectral indices n^. 
We have /(oo) ~ 200 for n</, = —1 {up — 1) model. This large value indicates that the 



^For the top-hat filter we have similar result /(oo) = 2.6. 



-19- 



PDF would become close to the Gaussian profile. The factor keeps f{R/r) < 3 for = 
—2 (up = —1) model. White (1999) investigated this model using N-body simulations and 
found that the PDF is skewed for Gaussian smoothing spanning a factor 5 in scale. This 
behavior is consistent with our result. For = —2.4 model, the factor f{R/r) changes 
only ~ 40% in the range < i? < oo. These dependence on the spectral indices seems 
reasonable considering that as increases, number of statistically independent region would 
also increases for a given smoothing volume and the PDF would become a more Gaussian 
like profile. 
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Fig. 1. — Nonlinear power spectrum A.{k) = PNL{k)k^ /{2t^'^) for n — —2 model (fig. la), 
n — —1.5 model (fig. lb), n — —1 model (fig.lc). In the right panels we plot the two-halo 
and one-halo terms separately. The total values (summations of two terms) are presented in 
the left panels. Dash-dotted lines represent the linear spectra. Thin solid lines correspond 
to the G-model, dotted-lines to the C-model, short-dashed lines to the LNl-model and long- 
dashed lines to the LN2-model. The thick-solid lines are obtained from the fitting formula of 
Peacock & Dodds (1996) given for Gaussian fluctuations. In the left panels the long-dashed 
lines (LN-2 model) and the dash-dotted lines (linear spectra) are nearly overlapped. 
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Fig. 2. — The mass functions and the bias parameters for various PDFs. Correspondence of 
hnes and models is same as figure 1. There are no halos with m > in the LN2-modeL We 
plot the mass function N{m) in the form N[m)m^, as the integral J N {m)m^ dm is important 
for our analysis. 




Fig. 3. — Nonlinear spectrum for non-Gaussian initial fluctuations with CDM spectra. We 
characterize the non-Gaussianity by the skewness parameter S. The ratio Pnl/Pl is given 
at /c = k]siL (fig 3. a) and k = O.GG/cat/, (fig 3.b). The solid lines represent the results from full 
dark-matter-halo approach and the dotted lines are obtained by the simple approximation 
described in the main text. 



-25- 




Fig. 4. — The scaling factor f{R/r) (see eq.[All]) as a function of Gaussian smoothing 
radius R. We have f{R/r) ~ 1 for R/r > 1 and f{R/r) ~ const for R/r > 20. 



